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ABSTRACT 






An experimental and analytical investigation was conducted to develop a calibrated 
torsional model of a three-cylinder, two-stroke diesel engine. A Detroit Diesel 3-53 engine 
was instrumented for time resolved measurement of cylinder firing pressures and high 
resolution near instantaneous shaft speed using a 720 and a 3,600 count per revolution optical 
encoder. Data were taken for three speeds and three torques for a total of nine conditions. A 
six degree-of-ffeedom torsional vibration model of the crankshaft, connecting rods, and 
pistons was developed. The nonlinear inertias, due to the reciprocating pistons, were included 
along with linear stiffness and damping. The equations of motion were numerically integrated 
over a cycle to obtain predicted response. The predicted response was compared to the 
measured response at the free end of the crankshaft. 
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I. INTRODUCTION 



A. MOTIVATION 

Diesel engines are widely used in both military and civilian applications. In the 
Navy, they are utilized for main propulsion and electric power generation. The reliable 
detection of faults both non-obtrusively and inexpensively could shift the focus of 
maintenance from regularly scheduled to condition based. This maintenance strategy has 
the potential to save costly, unneeded maintenance and can increase availability of diesel 
engines. 

Individual cylinder pressures of internal combustion engines are excellent 
predictors of many types of engine faults. Cylinder pressures are directly related to the 
torque output capability of an engine. Therefore, the contribution of torque from an 
individual cylinder is directly related to the maximum pressure that cylinder develops. 
Cylinders experiencing loss of pressure can be caused by faulty fuel injectors, worn 
cylinder lining, faulty piston rings, leaking valves, and blown or leaking cylinder head 
gaskets. 

There is a problem with measuring cylinder firing pressures in the field. While high 
quality pressure transducers can accurately measure cylinder pressures, they have a limited 
lifetime in the combustion chamber environment and are expensive. Many alternative 
indirect methods of cylinder pressure monitoring such as the use of strain bolts are being 
explored. 
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Another indirect method is to measure the variations in shaft speed. With this 
information and equations of motion describing the dynamics of the engine, the firing 
pressures are solved for numerically. In order to accurately predict cylinder firing 
pressures with shaft speed variation, an accurate torsional engine model representative of 
the physical system is required. 

B. STATE OF THE ART 

There have been numerous research efforts focused on extracting cylinder pressure 
information from crankshaft speed fluctuations. It is well known that these speed 
fluctuations are a direct result of the cyclic nature of the crankshaft torque waveform. 
This torque waveform is a direct result from gas pressure forces, due to combustion, and 
inertia forces, due to internal reciprocating masses. The challenge has been taking these 
measured speed fluctuations and passing them through a robust engine model to predict 
cylinder pressures. This literature review focused in areas of non-obtrusive cylinder 
pressure and torque prediction with emphasis on torsional vibration modeling. 

Wilson [Ref. 1] developed a torsional vibration model of a crankshaft with 
reciprocating components. With his methods, the individual inertia, stiffness, and damping 
between cylinder crankthrows and the rest of the shaft are readily determined. His 
methods are currently used by industry to develop torsional vibration models of 
crankshafts. Kabele [Ref. 2] developed a torsional vibration model of a diesel crankshaft 
which included piston and ring friction, connecting rod plane motion, hysteretic damping 
in the crankshaft, and cylinder pressure as a forcing function. The model was derived 



2 



from Newtonian Mechanics and the governing differential equations were non-linear 
making a numerical solution necessary. Once his model was developed, he compared 
measured and model predicted vibration amplitudes for various harmonic orders and 
concluded that correlation between the two was excellent. The measured vibration 
amplitudes were from a V-8 diesel engine. 

Citron et al, [Ref 3] did a computer simulation and produced cylinder pressure 
torque waveforms by using an elastic model of an engine, drivetrain system, and measured 
speed fluctuations. These speed fluctuation data were passed back through the model to 
determine the fluctuation waveform of both the total engine torque being developed and 
the cylinder pressure waveform which gave rise to it. 

Rizzoni [Ref. 4] proposed passing crankshaft speed fluctuations through an 
equivalent electronic circuit representing engine dynamics to predict cylinder misfires. He 
conducted experimental work on production spark ignited vehicles and confirmed that it is 
possible to apply this method with minimal hardware and computational overhead. 

Sobel et al, [Ref. 5] described a technique for measuring instantaneous crankshaft 
torque utilizing a non-contact ferromagnetic material to measure torsional stress in the 
crankshaft. The instantaneous torque of an internal combustion, four-cylinder engine was 
measured at the flywheel and compared with pressure signals during bench tests. 
Excellent correlation between this non-contact measured torque and values of the mean 
effective pressure for each cylinder were obtained. 

Brown and Neill [Ref 6] presented a non-contact method for determining 
simultaneously the pressure in each cylinder by employing a pattern recognition technique 
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to compare crankshaft speed fluctuations to reference patterns in a knowledge base. The 
non-contact method utilized was an interval timer and a magnetic sensor which timed the 
flywheel gear teeth as they passed. The experiment was carried out on a Detroit Diesel 
6V-92TA engine. The experimental results show that the method evaluates the cylinder 
pressures with an RMS error of less than six percent. But this requires extensive data 
bases for each condition. 

Bell [Ref. 7] concluded, utilizing the same Detroit Diesel 3-53 engine used in this 
research, that the variation of shaft speed held information which could be used to predict 
cylinder firing pressures. He also recommended the development of a torsional vibration 
model of the engine and flywheel system along with a more rigid optical encoder mounting 
for more accurate time resolved measurement. 

C. METHODOLOGY AND OBJECTIVES 

The engine used for this research is a Detroit Diesel Series 53 engine, model 5033- 
5001N. It is a three-cylinder, two-stroke engine. During one revolution, the crankshaft 
speed fluctuates as a result of cylinder pressure variations. These crankshaft speed 
fluctuations can be measured using an optical encoder. A torsional vibration model with 
nonlinear inertias, due to the reciprocating pistons, and linear stiffness and damping will be 
developed to represent the dynamics of the engine. The torsional vibration model will be 
calibrated by ensuring, as closely as possible, predicted response matches the measured 
response. Data from these measured speed fluctuations can then later be passed through 
the calibrated torsional model to determine cylinder firing pressures 
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The objectives of this thesis are to: 

1 . Instrument a diesel engine with a rigid optical encoder mount such that reliable time 
resolved measurements of crankshaft speed relative to the engine block can be obtained. 

2. Develop a torsional vibration model of the crankshaft including components connected 
to it (i.e. connecting rods, pistons, flywheels, and external loads). 

3. Operate the diesel engine at several speeds and applied torques to establish baseline 
data. 

4. Compare the model predicted response to the experimental measured response. 

D. ORGANIZATION 

Chapter II describes the development of the crankshaft torsional vibration model. 
Equations of motion are derived for a six degree-of-ffeedom torsional model and 
numerical solution procedure is presented. 

Chapter III describes the engine, dynomometer, instrumentation, and data 
acquisition procedures. 

Chapter IV presents results of data acquisition, the measured response, and the 
model predicted response. A comparison of measured to predicted response is made 
based on the results. 

Chapter V contains summary, conclusions, and recommendations for future work. 
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II. TORSIONAL VIBRATION MODEL 



A. PHYSICAL SYSTEM 

The engine used for this research is representative of diesel engines It is a three- 
cylinder, two-stroke engine. Engine operating conditions and design specifications are 
discussed in more detail in Section III. A. 

Figure 2.1 illustrates the engine layout with the crankshaft, flywheel, pistons and 
connecting rods. The front of the engine is to the left. The cylinders are numbered from 
front (1) to back (3). Connected to the power take-off shaft on the right is the load, 
which in this case is supplied by a water brake dynamometer. The firing order is 1-3-2 and 
the crankshaft rotates clockwise as seen from the front 




Figure 2. 1 Detroit Diesel 3-53 with crankshaft and major reciprocating 
components. The front of the engine is the left hand side. 
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The crankshaft crankthrows are the locations where the connecting rods are attached 
to the crankshaft. Crankthrows number one, two, and three are displaced 120 degrees. 
After the number one piston reaches top dead center (TDC), its highest point of travel, 
120 degrees of crankshaft rotation later the number three piston arrives at its TDC. 
Number two piston’s TDC is 240 degrees after number one. The power into the 
crankshaft is delivered at each crankthrow. 

The forces due to gas pressure in each cylinder can be taken as a force ‘-F’ on the 
piston and an equal and opposite force ‘F’ on the cylinder head. The force ‘-F’ is 
transmitted through the connecting rod to the crankthrow. The resultant force on the 
crankshaft multiplied by the eccentricity ‘R’ of the crankthrow creates the gas torque. If 
there were no friction, the gas torque would equal the shaft applied torque. Figure 2.2 
shows the relationship between gas pressure force and torque. 




where: F = cylinder pressure * area of the piston (lbf) 

R = crankthrow eccentricity (in) 

0 = angular position of crankshaft relative to TDC (rad) 
L = length of connecting rod (in) 



Figure 2.2 Relationship between gas pressure forces and gas torque into the crankshaft. 
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Gas torque into the crankshaft as a result of cylinder gas pressures is given by the 



following relationship: 



T i cy i=P icyiRsin( 6 )Ap 
T 2 cyl=P 2 cylRsin( 6 - 47t/3)Ap 
T 3c yi = P 3c y iRsin( 0 -27t/3)Ap 



( 2 . 1 ) 



where: P# cy i — Measured cylinder pressure (psia) 

R = Crankshaft crankthrow eccentricity (in) 
6 = Angular position of #1 crankthrow 
after TDC (rad) 

Ap = Cross-sectional area of the pistons (in 2 ) 



The torque applied to the crankshaft is somewhat less than the gas torques due 
mainly to friction on the cylinder wall (rings). The major components considered in the 
development of a torsional vibration model are the crankshaft, pistons, connecting rods, 
and flywheel. This physical system to be modeled is shown in Figure 2.3. 
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Figure 2.3 Crankshaft with a single reciprocating component. 



9 



1. Model Development 

The reciprocating components of the engine are very complex. However, to 
facilitate calculation it is useful to represent the complex system by a reduced order 
system. For this research, six rigid masses connected by sections of shafting with torsional 
stiffness and damping shall be used. This reduced system should retain the dominant 
dynamic characteristics of the original physical system that are of interest. The six degree- 
of-freedom torsional vibration model that is used to represent the physical system is shown 



in Figure 2.4. 
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Figure 2.4 Torsional vibration model of crankshaft and reciprocating components. 



In Figure 2.4, 0i represents the angular position of the crankshaft where the optical 
encoder attaches, 02 through 0 4 represents the angular position of crankthrows number 
one through three, 0 5 represents the angular position of the flywheel, and 06 represents the 
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